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Abstract 



Experiments on the motion of a particle on an inclined rough plane have 
yielded some surprising results. For example, it was found that the frictional 
force acting on the ball is viscous, i.e. proportional to the velocity rather 
than the expected square of the velocity. It was also found that, for a given 
inclination of the plane, the velocity of the ball scales as a power of its radius. 
We present here a one dimensional stochastic model based on the microscopic 
equations of motion of the ball, which exhibits the same behaviour as the 
experiments. This model yields a mechanism for the origins of the viscous 
friction force and the scaling of the velocity with the radius. It also reproduces 
other aspects of the phase diagram of the motion which we will discuss. 



I. INTRODUCTION 



A mixture of particles with different sizes can segregate under a wide range of conditions 
and due to various mechanisms (usually a combination) such as the "Brazil nut effect" 
shear 0, percolation 0, convection [Q, or surface flow ||^. Clearly the problems of 
mixing and segregation have important consequences for many industries. To understand 
this complicated phenomenon, it is easiest to study the effects of the various mechanisms 
separately, if possible. 

In this paper we are interested in segregation due to flow which can be seen, for example, 
in sand piles. When grains with various sizes flow on the surface, one observes that the largest 
grains find their way to the bottom of the pile while the smaller ones are stopped farther 
uphill. This segregation is caused by the roughness of the underlying surface on which the 
grains are flowing. In a much simplified picture, one can consider the bulk of the sand pile 
as providing a rough substrate on which the surface grains are flowing and segregating. For 
experimental purposes, one can carry the simplification farther and consider an inclined 
plane made rough by sticking beads of a given radius, r, for example glass or sifted sand, on 
contact paper and attaching it to the plane. This was done in a series of experiments @-|^ 
where various properties of the motion of balls down this plane were studied. We refer the 
reader to the references for details of the experiments and results. Here we will review some 
of the results before proceeding to our theoretical model. 

Three regimes were found for the motion of the ball down the plane: A) A sticking 
or pinning regime where the ball comes to a sudden stop after travelling some distance, B) 
a regime where the ball attains a steady state with a constant (on average) velocity inde- 
pendent of the initial release velocity, and C) a jumping regime where the ball experiences 
big bounces and does not achieve a steady state on the 2 meter long plane used in the 
experiments. One expects that for a long enough plane, the ball will always reach a steady 
state, but even so, the nature of the motion in this regime appears to be different from 
that of region B. However, because of experimental difficulties, this regime remains mostly 
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unstudied. 

In region B, it was found that the constant (average) velocity satisfies (v) oc i?"sin(^) 
||1CI|| . This is a surprising resuh because it means that the force of friction acting on the ball 
is viscous, i.e. proportional to the velocity l§-p|. A straightforward argument suggests that 



the friction force should be quadratic in the velocity leading to a ^sm{6) behaviour: It is 
clear that the deceleration due to collisions is proportional to the number of collisions per 
second times the velocity loss per collision. The velocity loss per collision is proportional 
to the velocity itself, and one can argue that the average number of collisions/second is the 
velocity of the ball divided by the average distance between surface beads. This gives a 
friction force which is quadratic in the velocity, contrary to what is seen experimentally. 
Two assumptions that would lead to viscous friction are 1) the velocity loss per collision 
is a constant independent of the velocity, or 2) the number of collisions per second is a 
constant independent of the velocity. However, neither of these assumptions can be justified 
physically and the explanation should be sought elsewhere. 

It is tempting to draw an analogy between the motion of a ball down a rough plane with 
peaks and troughs and the motion of particles in random potentials using a Langevin 
equation to describe the dynamics. However, such an approach cannot explain the viscous 
dissipation because the viscous friction term is put in as an assumption from the start. In 
addition, even if one does assume the viscous force and pursues this approach, the scaling 
of the velocity with R, i.e. with the roughness which in this approach characterizes the 
random potential, is not correctly reproduced. The missing ingredient will be discussed 
in section III. A different approach was taken in |^ where a Langevin equation, with the 
viscous dissipation and the velocity scaling with R built in, was used to study other aspects 
of the motion such as the stopping distance and its dependence on various parameters. 

The approach we will take here is different still. We will consider the one dimensional 
case, i.e. discs of radius r are randomly stuck on a line inclined at an angle 6 with the 
horizontal, and on which is released a disc of radius R. We then consider the equations of 
motion of the disc, R, as it moves down the line, and make as many simplifying assumptions 



as is reasonable, eventually turning the deterministic equations into stochastic ones. Our 
goal in this is not to have a detailed microscopic agreement with experiments: That would be 
an unrealistic hope with a very simple model and is more the domain of detailed numerical 
simulations. Rather, our goal is to have a bare-bones stochastic model that agrees with 
experiments and reproduces various properties of the motion such as the viscous friction 
and the scaling of the velocity with R. With this model, we succeeded in finding a simple 
mechanism for the viscous force and the various scaling laws. In all this geometry plays a 
crucial role. 

In section II we will construct the model clearly stating all our assumptions, and in 
section III we will compare our results with the experiments and make some comments. 
Section IV contains our conclusions. 

II. THE MODEL 

Figure 1 shows the geometry of the system. Discs of radius r are put in a straight line, the 
distance between the surfaces of two neighboring discs being 2e, where e is a random number 
between and a maximum value, cm- This line is inclined an angle 9 with respect to the 
horizontal, and we will take the direction parallel (perpendicular) to it as the x (y) direction. 
Therefore, the gravitational acceleration in the x {y) direction is gisin(^) (— gfcos(6')). The 
big disc, radius i?, moves from right to left with a velocity v and x and y components Vx and 
Vy, making an angle with the line of discs, sin(0) = As shown in Fig. 1, 7 is the 

angle between the line perpendicular to the line of discs, and the line connecting the centers 
of the big and small discs at the point of contact. We will call this the angle of contact 
which, it should be emphasized, is different from the angle of incidence. Our sign convention 
is that when the collision is on that side of the small disc which faces the approaching big 
disc (as in the figure) 7 < 0, and when it is on the other side, 7 > 0. 

Our first assumption is that we can ignore the rotation of the big disc as it moves down 
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the line. Clearly this assumption is not justified if we are aiming at detailed microscopic 
comparison with experiments. For example, neglecting rotation means that Vt-, the tangential 
velocity of the point of contact during a collision is due entirely to the translational velocity of 
the disc. This is important quantitatively because ^t-, the tangential coefficient of restitution, 
depends on Vt (see below). However, molecular dynamics (MD) simulations have shown 
|l^ that if we prevent the disc from rotating as it bounces down the line and have only 



dynamic friction acting as the tangential force, the qualitative features of the motion, like 
the 9 dependence of the velocity, do not change. Since we are interested here in the scaling 
properties of the motion, we will neglect the rotation. Therefore, with this assumption, the 
only effect of the size of the disc is to contribute geometrical constraints as we will see below. 

Our second assumption is that after a collision, the velocity components normal and 
tangential to the small disc at the point of impact are related to the corresponding velocities 
before the collision by 

Wn\=l^n\Vn\, (lo) 

v[ = HtVt, (lb) 

where /i„(/it) is the normal (tangential) coefficient of restitution. There are absolute value 
signs in eq.(la) but not eq.(lb) because the coefficient of normal restitution is always positive 
while that for tangential restitution can be negative (see eq.(9a,b)). The angle of incidence, 

(not shown in the figure), is of course given by cos(aj) = Ifnl/liJI, and sin(aj) = 
while the angle of reflection, a^, is given by cos(ar) = '|, and sin(ar) = '|, 

where the primes denote values just after the collision. Now consider the big disc going 
down the inclined line colliding with the small discs, and let us examine the kth collision. 
Let the x and y velocities just before (after) the kth collision be Vx{k) (v'^ik)) and Vy{k) 
{Vy{k)). It is easy to see from these deflnitions and the simple plane geometry that 

Vn{k) = \vy{k)\cos{'j) - Vx{k)sm{'y), (2a) 

vt{k) = |fy(A;)|sin(7) + Vx{k)cos{-f), (26) 
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keeping in mind our convention for the sign of 7. The x velocity after the kth colhsion is 
v'^{k) = \v '|sin(a.r + 7) which, after expanding the sine and using the above definitions of 
the angle of reflection, gives 

v'^{k) = HtVtCOsi'j) + HnVnSm{'j). (3) 

Combining this with eqs.(2a,b) allows us to express the x velocity just after the kth collision 
in terms of the x and y velocities just before the collision as 

v'^{k) = v^{k)(^fj,tCos^{-f) -/x„sin^(7)^ + |fj^(A;)| (^/i* + /in^sin(7)cos(7). (4a) 

A similar argument for Vy{k) gives 

v'y{k) = -w^(A;)(^/ii + /i„^sin(7)cos(7) + |wy(A;)| (^yU„cos^(7) -/iisin^(7)y (46) 

Our next assumption, which is supported by MD simulations [|T^,|T^, is that the motion 
of the disc is composed mainly of a series of small bounces with very little or no rolling which 
we will, therefore, ignore. This means that the kth collision will cause the disc to bounce 
and spend a time (5r(/c) in the air during which it will experience x acceleration gsm{6). 
Therefore, its x velocity just before the (A; + l)th collision is Vx{k + 1) = Sr^k) gsm{9) + v'^{k) , 
which, combined with eq(4a), gives 

v^{k + l) = 6T{k)gsm{e) +v^{k)( ^u^cos^ (7) -/i„sin^ (7) ) + l^y (^) I ( l^t+l^n ) sin(7)cos(7) . (5a) 



What happens to the perpendicular velocity, Vy7 If the disc, after the collision shown in fig. 
1, bounces up and lands on another small disc but at the same y value, its Vy just before the 
new collision is identical to that just after the previous collision. In general, however, the y 
value will be different because the big disc will land on various different points of the small 
discs. However, since we are typically interested in R/r > 3, we see that the variations in 
y are very small. We can, therefore, make the simplifying assumption that Vy{k + 1), the 
y velocity just before colhsion {k + 1) equals Vy{k), the y velocity just after collision k |]14 . 
We then have 

Vy{k + 1) = -Vx{k)(^nt + /i„^sin(7)cos(7) + |t;3/(A;)| (^/i„cos^(7) - fxtsin^ (l)^ ■ (56) 



So, eqs.(5a,b) express the x and y velocities just before the {k + l)th coUision in terms of 
those just before the kth coUision. Although not written explicitly, to simplify the notation, 
it should not be forgotten that the contact angle 7 appearing on the right hand side of these 
equations is actually "y{k), the angle for the kth collision. These equations are deterministic 
since, given the initial conditions, one can, in principle, calculate for all subsequent collisions 
the two terms that have not yet been specified, ST{k) and j{k). Such a detailed microscopic 
approach is not our goal in this paper. Instead, we want to transform eqs.(5a,b) into stochas- 
tic processes by making 7 a stochastic variable. To do that we need its distribution which 
must be based on the underlying microscopies of the collisions. To simplify the calculation 
of this distribution, we assume that the disc, R, collides with a given small disc only once. 
It is clear from fig. 1 that, for a given incident velocity, v, or equivalently, for a given angle 
4>, the largest (positive) contact angle, jmax{k), is obtained when the big disc collides with 
the small disc while tangent to the hne L2. For this case, ^max{k) — (j){k). On the other 
hand, the smallest {i.e. the largest negative) 7mm(^) is obtained for the case shown in the 
figure, with the big disc tangent to line Li. Notice that Li is tangent to the disc just to 
the right of the disc involved in the coUision. This shows how discs cast velocity dependent 
"shadows" on their neighbors thus restricting the area available for a collision with a given 
bead. It is straightforward to show that 

sm{jrnin{k)) = sm{4>{k))cos{a{k)) — cos{4>{k))sm{a{k)) , (6a) 

where 

cos(a(A;)) = 1 - 2^^^sin(0(A;)), (66) 
r + R 

and where 2(r + e) is the distance between the center of the disc under collision and that to 
its immediate right. Recall that e is a uniformly distributed random number between and 
some maximum value em- This gives the range of the stochastic variable, j{k), to be between 
a maximum value, (f){k) (where sin(0(A;)) = Vy{k)/\v{k)\), and a minimum value given by 
eqs(6a,b). But we still need its distribution. To this end, we assume that for a given collision, 
the point of disc R that is closest to line Li (in fig. 1 it actually touches it) is equally likely 
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to be anywhere between Li and L2. The position of this point determines the value of the 
contact angle 7, and therefore knowing its probability distribution (which we assumed to be 
uniform) gives the probability distribution for 7. Using this, it is straightforward to show 
from the geometry that 7 is given by 



where 

r 



y — uniform random number e 



-2(r + e),- 



(76) 



.sin(^) 'sin(^). 
One sees that the distribution for 7 is a function of the geometry {R, r, and e) and 

the incoming velocities {vx and Vy) which determine how much of the surface disc area is 

available for collision. The distribution given by eqs. (7a,b) is shown in fig. 2. 

Equations (5a,b) and (7a,b) define most of the stochastic model. We still need to specify 
ST{k), the time between collisions, when energy is fed into the system. Perhaps a temptingly 
simple hypothesis is that 5T{k) — d/v'^{k) where d is the average distance between the centers 
of the line discs. However, as explained in section I, it is easy to see that this gives a friction 
force proportional to v^, and not Vx- An alternative argument is to say that for a perfectly 
smooth surface, the time between bounces is given only by Vy {ST{k) — 2vy{k) / gcos{9)) , and 
assume the same holds for the rough plane. Integrating the stochastic equations of motion 
with this assumption produces two types of solution depending on the parameters used: 
Either the ball energy is continuously dissipated until it comes to a complete stop, or the 
ball makes bigger and bigger jumps accelerating all the time and never reaching a stationary 
state. Clearly neither of these two possibilities for ST{k) is realistic. Since the surface is 
rough, both Vx and Vy play a role in determining 5T{k). As we shall see below, the balance 
between Vx and Vy plays a crucial role in determining the properties of motion. 

To motivate our choice for ST{k) note, from eqs. (4a,b), that a collision with a negative 7 
transfers velocity from the x to the y direction, while the opposite happens for a positive 7. 
Furthermore, for a negative 7 collision, the big disc has to jump over the small disc in order 
to continue its motion, while for a positive 7 collision the big disc can always continue its 
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motion. For these reasons, we assume that for a coUision with 7 negative (positive), ST{k) 
is determined only by Vy {v^). In other words, for ^{k) < 



and for ■j{k) > 



Sr(k) = (86) 



Note that in eq. (8b) we are assuming that the distance travelled to the next collision is well 
approximated (on average) by 2(r + e). Of course in reality, v^; and Vy together determine the 
time between collisions. This can be calculated but would not reveal the underlying reasons 
for the viscous force which we are trying to explain. Equations (8a,b) give a simplification 
that attempts to understand the separate roles of the parallel and perpendicular velocities. 

The last ingredient in our model is the experimental observation that the tangential 
coefficient of restitution is strongly dependent on the angle of incidence |T5|, given by the 



normal (tangential) components of the velocity, Vn (ft) (see eqs.(2a,b)). It was found in 
reference |15| that the tangential coefficient of restitution is well parameterized by 

V (k) 

/it(A;) = 1-3.5(1 +/i„)/i^^, (9a) 

vt{k) 

in collisions that involve gross sliding, whereas in collisions that do not 

Mk) = (9&) 

In these equations, fi is the coefficient of friction, and f3 is the tangential coefficient of 
restitution in the absence of sliding. The normal and tangential velocities are given by 
eqs.(2a,b). We mention here that one consequence of our ignoring the rotation of the big 
disc is that the tangential velocity is given only in terms of the translational velocities and 
Vy. However, rotation modifies the tangential (but not the normal) velocity and therefore 
the effective fit for each collision. This effective change in fit will lead to a different average 
velocity, but not to a qualitatively different behaviour. This is one justification for ignoring 
rotation. 



Now, our stochastic process is completely defined by eqs. (5a,b), (7a,b), (8a,b) and 
(9a,b). This stochastic process is nonlinear and cannot be solved exactly. However the 
equations can be easily iterated numerically, and the statistical properties of the motion of 
the disc down the line studied in detail. The initial conditions are random choices for Vx{l), 
Vy{l) and 7(1) from which everything else follows simply by generating the random number 
y. If a collision results in a negative value for Vx we consider the disc to have been pinned. 

The equations defining our model can be further simplified by noting that all the angles 
are very small allowing us to expand all trigonometric functions to second order in 7. How- 
ever, the resulting equations will still be nonlinear and require a numerical integration. The 
utility of such an expansion would be for approximate solutions, e.g. mean field. 

III. RESULTS 

In this section we compare the results of our model with the experiments. The first test 
of the model is whether it can reproduce the three experimentally found phases the 
pinning phase A, the constant velocity phase, B, where (vx) oc sin(^), and the "jumping" 
phase, C. Figure (3a) shows (vx) as a function of sin(6') for three values of R, the radius of 
the big disc. We clearly see regions that are linear in sin(^) through which we have fitted 
straight lines. Although we will not attempt to match our results exactly with experimental 
values, it is beneficial to see if the theoretical results are in reasonable general agreement 
with experiments. The linear regions stretch from about 4.5 to 14.5 degrees (the precise 
values depend on R, and for a given R they depend on the coefficients of friction and normal 
restitution), and the values of (vx) are of the order of 5 to 15 cm/sec, all in very good general 
agreement with the experiments. Also in agreement with experiments, we found that the 
average velocity in region B does not depend on the initial velocity with which the ball was 
released. This is an important property of the motion in this region and does not apply in 
region C. When we scale the curves in fig. (3a) by a factor of R~^'^^ we get fig. (3b) which 
shows excellent scaling in the linear region, again in agreement with the experiments. The 
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exponent we found, —0.49 differs from the experimental value —1.4. Changing the 

parameters of the calculation, i.e. the values of the coefficients of restitution and friction, 
changes the exponents only a little and does not bring it close enough to the experimental 
value. We think that this difference could be due to the fact that the calculation is done 
in one dimension while the experiments are two dimensional. It might also be that our 
approximations, while maintaining the qualitative features of the physics, have changed this 
exponent. This insensitivity of the exponent to the values of the coefficients of restitution 



and friction explains the experimental result ||I6| that for a given value of 9, R, and r, the 



average velocity is practically independent of the material of the plane or the ball In 
other words, having fixed 9, R, and r, balls of plastic, glass and steel, whose coefficients 
of restitution and friction are not very different, were found to have the same terminal 
velocity. Furthermore, this terminal velocity did not change when the surface of the plane 
was changed from glass beads, radius r, to sand with average grain radius r. 

To check how linear these plots are, we show in fig. (4) (vx) —Vfu where vju is the value of 
the velocity given by the straight line fit to each of the curves. We see that for sin(6') between 
about 0.1 and 0.25 the deviation is extremely small and the linear fit is excellent. For sin(^) 
larger than about 0.25 we see that the (vx) increases faster than linearly. We take this as 
the definition of region C, which for brevity we will call the "jumping" regime because, as 
mentioned in the introduction, the motion is dominated by large jumps as compared with 
the small bounces in region B. We see from fig. (4) that the largest disc, dashed line, enters 
the jumping regime earlier and faster than the smaller ones. This is in agreement with 



experimental phase diagrams and with numerical simulations [|18],|T2|. Experiments 



have also shown that as the ball aproaches the jumping regime, C, it exhibits intermittent 
behaviour where it still has a constant average velocity but it shows bursts of acceleration 
and deceleration. Figure (5) shows a plot of distance travelled as a function of time for 
sin(6') = 0.53 from our model. Bursts of acceleration and deceleration are clearly seen. This 
intermittency gets much more dramatic as one gets closer to the jumping regime (increasing 
9) before the motion gets completely destabilized. The origin of a burst of acceleration is 
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a particularly big bounce, due to high speed and a large negative 7, where the disc spends 
a long time in the air, gaining a lot of energy which it has difficulty losing. Then comes 
another big shock that dissipates a lot of the energy, and the disc regains it former slower 
velocity. This mechanism can be inferred from the sharpness of slope changes [i.e. velocity 
changes) in fig. (5). 

Figure (4) also shows that as 6 is decreased, the average velocity decreases and eventually 
exhibits a sublinear dependence before it enters the pinning regime. To study the transition 
between the pinning regime. A, and the constant velocity regime, B, we release several discs 
for each value of 6 and R (keeping r fixed) one after the other and with slightly different 
initial velocities. We consider the system in the pinning regime when at least half the 
discs are stopped. As is clear from fig. (6), the pinning regime quickly shrinks {i.e. moves 
towards smaller 6), and then disappears as R is increased. This agrees with the experiments, 
although in our calculation region A disappears faster. 

Experiments have indicated that (vx) does not depend on the ratio R/r alone. 

For example, Jan et al ||TP[ showed that changing R while keeping the ratio R/r fixed (they 



kept R = r), {v^) grows like ^/R. Figure (7) shows {vx)R~'^'^ as a function of sin(6') for 
R = 3,4.5,6,7.5,9mm and R/r fixed at 6. The figure exhibits excellent data collapse with 
an exponent that agrees with the experimental value. The same scaling is observed in MD 



simulations [12 



So, we see that the simplified stochastic model presented here has reproduced remarkably 
well all the features of the experimental results. But what is the mechanism behind the linear 
dependence on sin(^) and the scaling with R7 As mentioned in section II, if one argues that 
the time between collisions, 6t, is set by the parallel velocity, v^, one gets a very stable 



dependence of the velocity on ysm{6) implying a friction force that is proportional to the 
square of the velocity. The faster the disc moves, the shorter the time interval during which 
energy is fed into the system. If, on the other hand, one assumes, as in ballistic flight, that 
6t is determined only by the perpendicular velocity, Vy, then one gets very unstable motion. 
The faster the disc, the longer the time interval during which energy is fed into the system 
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and the more difficult it is to dissipate it. This drives the instabihty. In our simple model, 
the time of ffight is given by both and Vy depending on the side of the line beads that is 
impacted. 



fig. (8a) for sin(^) = 0.03. This gives a dependence on sin(^) that is very close to a square 
root, as is seen for the same angle in fig. (3a,b). The reason for this is that at such small 
angles, Vy is so small that the motion is mostly parallel to the line, i.e. in the x direction. 



which as mentioned above, gives a Ysin(^) dependence. For large values of 6, the motion is 
very bumpy, Vy is quite large and it dominates the contribution to 5r, as is seen in fig. (8c) 
for sin(^) = 0.4. At this value, we see that {v^) increases faster than linearly with sin(^), 
fig. (3a,b), as the system is headed for the jumping regime, C. In between the two regions 
of sublinear and superlinear dependence on sin(^), there is a region of competition between 
the Vx and Vy contributions to 5t as is seen in fig. (8b) for sin(0) = 0.2, which is in the heart 
of the linear regime. It is this competition between the destabilizing influence of Vy and the 
strong stabilizing influence of that gives rise to the effective viscous friction and linear 
dependence on sin(6'). 

What about the scaling observed in figs. (3b) and (7)? We found that if we use a 
tangential coefficient of restitution, /it, that is constant instead of as in eq. (9a,b), the 
scaling exponent for fig. (3a) is 0, i.e. all the curves always overlap completely in the 
linear regime. Furthermore, we found that with this assumption, it is the smaller discs that 
enter the jumping regime, C, ffist, in clear contradiction with the experiments and intuition. 
Therefore the experimentally based eqs. (9a,b) are crucial to obtaining the correct behaviour. 
On the other hand, we found that taking /if to be constant does not change the behaviour 
in fig. (7) {i.e. when R/r =constant) by much: the exponent is changed from 0.5 to 0.43 
and the scaling is still clearly visible. This interesting result indicates that although, in our 
model, the exponents for figs. (3b) and (7) are very similar, the two scaling effects appear 
to have different origins! 

It is worth emphasizing that geometry is playing a dominant role in this model. The 



For small inclination angles, our model gives 5t that is dominated by Vx as shown in 
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distribution of the stochastic variable, 7, given by eqs. (7a,b), is dominated by the velocity 
dependent shadows cast by discs on their neighbors, constraining the cross section available 
for collisions. In addition, this distribution is clearly crucial in deciding which of eqs. (8a) 
or (8b) determines the time of flight whose role is discussed above. Because of all this, we 
see that the allowed collisions and the velocities are strongly dependent on the geometry 
which also means that, because of eq. (9a), is also dependent on the geometry. 



IV. CONCLUSIONS 

We have presented in this paper a stochastic model for the motion of a disc down an 
inclined line of smaller discs, separated by random distances. The stochastic equations are 
based on the original deterministic equations of motion to which we have added some simpli- 
fying assumptions. Our two main assumptions are that we can neglect rotation (supported 
by MD [1^), and that the motion of the disc is a series of small bounces, with one bounce 
per line disc (supported by recent experimental results on the two dimensional plane and 



by MD simulations [Q). This simple model, where geometry plays a very important role, 
accurately reproduces the features of the experimental data, and explains the origin of the 
viscous friction force. This is seen to be the result of competition between ballistic motion 
(with large Vy determining the the time of flight) and motion parallel to the plane (where 
Vx determines the time of flight). We also emphasized the importance of the dependence of 
the tangential coefficient of restitution on the angle of incidence. 

The next phase is to compare the results from this model with experiments currently 
being done on the stopping distance of the balls in regions A and B. The properties of the 
stopping distance in region B have very important industrial applications in the segregation 
of grains with different sizes. We will also compare with experiments the results of our model 
for the longitudinal dispersion of the velocity and its dependence on the sizes of the balls 
and the inclination angle. For the longer term, we would like to generalize the model to a 
two dimensional plane, where we can make more quantitative comparisons with experiments 
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including results for the transverse motion which is absent on a one dimensional line. 
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FIGURE CAPTIONS 



FIG. 1 Shows the geometry of the inchned hne of discs, r, on which we release a disc R. The 
coUision depcited is that for the largest negative value of 7, the contact angle. 

FIG. 2 A histogram showing the distribution of the stochastic variable 7 (in radians) given 
by eqs. (7a,b). Recall that in our convention, 7 < if the collision is as shown in fig. 
1, and 7 > if it is on the other side of the perpendicular, sin(^) = 0.3, R — 0.3cm, 
r = 0.05cm, e = 0.01cm, = 0.9, = 0.142 and (3 = 0.45. See eqs. (5a,b), (7a,b), 
(8a,b) and (9a,b) for definitions of the parameters. 

FIG. 3a Shows the dependence of (vx) (in cm/sec) on sin(^^). There is a clear linear region 
in agreement with experiments. R — 3,4,5mm for □, A, y, /i = 0.142, = 0.9, 
P = 0.45, r = 0.5mm. 

FIG. 3b Same as fig. (3a) but with the vertical axis scaled by 

FIG. 4 Deviations of < > from the straight line fits of fig. (3a). R — 3,4,5mm for the 
dot-dash, sohd, and dashed lines. 

FIG. 5 Distance travelled (in meters) as a function of time (in seconds) for a disc of radius 
R = 4mm and an inclination angle sin(^^) = 0.53. All the other parameters are as in 
fig(3a,b). The slope gives the velocity parallel to the inclined line and exhibits bursts 
of acceleration and deceleration. 

FIG. 6 Fraction of pinned discs as a function of the inchnation angle. R — 2, 2.2, 2.4, 2.6, 
2.7mm for □,<0>, A,V;+- The other parametrs are as in fig (3a,b). For each data 
point 500 discs were released. 

FIG. 7 Scaling of {v^) with R when R/r is kept constant. The data collapse is for R = 
3, 4.5, 6, 7.5, 9mm and the same parameters as in fig(3a,b). 
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FIG. 8a Histogram of times between bounces. The shaded histogram shows 5t for coUisions 
at positive 7, i.e. where the time of flight is determined by v^. The other histogram 
shows the case for negative 7, i.e. St is determined by Vy. Sin(6') = 0.03, the other 
parameters are as in fig. (3a,b). 

FIG. 8b As in fig. (8a) but with sin(^) = 0.2. 

FIG. 8c As in fig. (8a) but with sin(^) = 0.4. 
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